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Se.c.tion  I 

PLASMA  THEORY 


A.  DRIFT  CYCLOTRON  INSTABILITY 
J.  K.  Lee 

We  use  the  dispersion  equation  solver,  ROOTS*,  to  obtain  the 
results  of  the  linear  theory  (growth  rates,  etc.)  as  a guideline  of  and 
a comparison  with  future  simulations  of  density-gradient-induced  instab- 
ilities, namely,  drift  cyclotron  (DC),  drift  cyclotron  loss  cone  (DCLC)  , 
or  lower  hybrid  drift  (LHD)  instabilities. 

For  the  drift  cylclotron  instability,  ROOTS  with  the  local 

approximation  provides  a dispersion  diagram  virtually  identical  to  Fig.  1 

2 

of  Gary  and  Sanderson's  local  Vlasov  theory.  In  addition  to  this, 

ROOTS  provides  nonlocal  solutions,  which  will  be  useful  in  most  simulations 

with  steep  density  gradients.  Some  details  will  appear  in  the  next  QPR. 

We  found  the  ROOTS  results  in  the  linear  local  dispersion 

solutions  to  agree  with  Cohen  eX  aX ROOTS  predicts  the  linear 

marginal  stability  at  L /p.=5.0  for  m./m  = 100,  id  /w  =1  (where 
1 dn  n 1 1 e pe  ce 

*"nE77dx*  's  t*ie  9u'd'n9  center  density,  p.  is  the  ion  gyroradius, 

to  and  to  are  the  plasma  and  cyclotron  frequencies  of  electrons), 
pe  ce 

Detailed  dispersion  diagrams  for  a choice  of  future  simula- 
tion parameters  will  be  given  in  the  next  QPR. 


1M.  J.  Gerver,  "ROOTS,  A Dispersion  Equation  Solver",  UCB/ERL  Memorandum 
No.  M77/27  (Oct.  31,  1976). 

2S.  P.  Gary  and  J.  J.  Sanderson,  "Density  Gradient  Drift  Instabilities: 

Oblique  Propagation  at  Zero  Beta",  Phys.  Fluids  2^,  pp.  1181-1187  (July,  1978). 

3B.  I.  Cohen,  N.  Maron  and  T.  D.  Ro^nlien,  "A  Local  Theory  of  Nonlinear  Ion 
Dynamics  in  a Drift-Cyclotron  Mode",  OA5,  Fusion  Theory  Conference  (197G) . 
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B.  PLASMAS  WITH  FIELD  REVERSAL* 
Douglas  Harned 


(1)  MAGNETIC  RECONNECTION 

The  problems  of  magnetic  reconnection  have  been  examined  to  see 
if  an  analogy  can  be  made  between  the  bow  shock  and  the  x-points  in  a field- 
reversed  mirror  configuration.  The  situations  appear  to  be  considerably 
different  because  of  the  requirement  for  particle  flow  across  magnetic 
field  lines  in  the  reconnection.  While  the  solar  wind  provides  such  a 
flow  in  the  magnetosphere,  the  relatively  low  particle  density  in  the  reqions 
surrounding  the  x-points  in  the  field-reversed  mirror  suggest  that  reconnec- 
tion will  not  occur  at  a significant  rate.  There  is  the  possibility,  however, 
that  in  a fusion  plasma,  alpha  particles  could  provide  a significant  source 
for  the  reconnection  process. 

Additionally,  some  facility  was  developed  with  the  particle 
2 

trajectory  code  TIBRO-X.  It  would  be  possible  with  this  code  to  investi- 
gate single  particle  behavior  in  simple  reconnection  schemes  by  replacing 
the  gravitational  force  with  an  electric  one.  New  models  similar  to  the 
current  sheet  models  of  Speiser^  could  be  examined. 

(2)  STABILITY  OF  REVFRSEV-FIELV  CONFIGURATIONS 

An  investigation  has  been  started  to  analyze  the  stability  of 

infinite  layer  (infinite  in  z)  reversed-f ield  configurations  to  modes 

having  k =0.  One  motivation  for  this  study  has  been  tht  recent  results 
z 45 

of  simulations  at  Livermore  and  Cornell  showing  the  instability  of  such 

*This  section  was  in  Section  II:  Simulation  in  the  past  QPR,  but  should  have 
been  in  Section  I:  Theory. 

Wasyliunas,  V.  M.,  "Theoretical  Models  of  Magnetic  Field  Line  Merging,  I",  Rev.  of 
Geophysics  and  Space  Physics  1J,  303-336  (1975). 

2J.D.  Hansen  and  R.  S.  Devoto,  "TIBRO-X:  A Code  to  Compute  Trajectories  of 
Charged  Particles  in  Magnetic  Fields",  UCRL-52189  (1976). 

3Speiser,  T.  W. , "Particle  Trajectories  in  a Model  Current  Sheet,  Based  on 
the  Open  Model  of  the  Magnetosphere,  with  Applications  to  Auroral  Particles", 

J.  of  Geophysical  Research  70,  1717-1728  (1965). 

4Byers,  J.  A.,  private  communication  (1978). 

5Friedman,  A.,  private  communication  (1978). 
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a cylinder.  We  are  beginning  by  analyzing  a slab  of  plasma  containing  a thin 
beam  (the  limit  of  an  infinite  radius  cylinder).  Hopefully  this  will  provide 
insight  into  the  stability  of  such  configurations  and  may  help  determine 
whether  the  instabilities  seen  elsewhere  are  real  or  numerical  in  nature. 


SzctLon  II 
SIMULATION 

A.  QUASIL INEAR  SIMULATION  OF  A MIRROR  MACHINE 
Dr.  Y.  Matsuda  with  Dr.  H.  L.  Berk  (LLL) 

A report  is  in  preparation  titled  "Quasi  I inear  Simulation  of 
the  Drift-Cyclotron  Loss  Cone  Instability". 

B.  SELF-HEATING  OF  A Id  THERMAL  PLASMA 
A.  Peiravi 

This  project  was  completed  this  past  quarter.  The  results 
were  reported  at  Monterey  (See  Section  V)  and  an  ERL  Report  is  completed, 
to  be  mailed  soon . 

Dr.  Hirotada  Abe  (Dept,  of  Electronics,  Kyoto  University)  has 
kindly  called  our  attention  to  his  work^  dealing  with  the  same  subject. 

We  recommend  this  work  to  our  readers. 

C.  ONE  DIMENSIONAL  ESI  CODE  FOR  SHEATHS 
Yu-Jiuan  Chen  (Prof.  C.  K.  Birdsall) 

No  special  work  to  report  this  quarter. 

D.  PARTICLE  TRAJECTORIES  IN  A CUSP  FIELD 
Yu-Jiuan  Chen  (Prof.  C.  K.  Birdsall) 

No  ap0«  ial  k t%»  I d|’.'l  I 1 M * •)»»*'  I 

1H.  Abe,  J.  Miyamoto  and  R.  Itatani,  "Grid  Effects  on  the  Plasma  Simulation 
by  the  Finite-Sized  Particle",  J.  Comput.  Phys.  Jj),  pp.  1 31*- 1 ^9  ( Oct .,  1975) 
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CODE  DEVELOPMENT  AND  MAINTENANCE 


r r 


A.  ESI  CODE 
Dr.  W.  M.  Fawley 

(1)  Adaptation!,  ofi  ESI  to  LI ie.  Po&t-Pnoaal>i>oK  "ZEP" 

I 

We  have  decided  to  adapt  the  much-used  ljd  electrostatic 
code  ESI  to  produce  time  history  files  that  can  be  read  and  manipulated 
by  the  post-processor  program  ZED.  ZED  is  normally  used  with  the  output 
of  the  2id  code  ZOHAR  and  permits  the  user  to  investigate  plasma  simu- 
lation time  histories  of  quantities  such  as  electric  potential  and  poten- 
tial particle  moments  through  means  such  as  Fourier  analysis  and  cross- 
correlation. It  is  necessary  to  change  ESI  to  run  under  the  CHATR  compiler. 

In  the  next  QPR,  we  plan  to  give  extensive  details  concerning  these  changes. 

(2)  A.  B.  Langdon  has  changed  the  loader  for  B/0.  Consult  ES1REP 
in  ESI  LIB. 

B.  EMI  CODE 
Dr.  W.  M.  Fawley 

ZED  will  also  be  set  up  for  EMI  in  the  near  future. 

C.  ELECTROSTATIC  2±d  CODE  DEVELOPMENT  (EZOHAR) 

W.  M.  Nevins,  Dr.  Y.  Matsuda 

This  work  is  nearly  complete,  as  reported  at  Monterey  (see 
Section  V)  and  is  being  prepared  for  publication,  with  the  title  "A  New 
Particle  Simulation  Model  for  a Bounded  Plasma",  by  Nevins  and  Matsuda. 

An  initial  application  of  the  code  will  be  given  in  a paper 
titled  "Simulation  of  the  Drift  Cyclotron  Instability"  by  Matsuda,  in 
prepa rat  ion. 


D.  ESI  + EFL  CODE 

J.  K.  Lee  (Prof.  C.  K.  Birdsall) 

A paper  titled  "Particle-Fluid  Hybrid  Codes  Applied  to  Beam- 
Plasma,  Ring-Plasma  Instabilities"  was  presented  at  Monterey  (see  Section 
V).  The  abstract  is  as  follows: 


PARTICLE- FLUID  HYBRID  CODES  APPLIED  TO  BEAM-PLASMA,  HINT, -PLASMA  IHSTABILITUS 

J.  K.  Lee  and  C.  K.  BlAdiaL'l 

Electronics  Research  Laboratory 
University  of  California,  Berkeley  CA  94720 

When  a beam-plasma  type  instability  is  simulated  using  a particle 
code,  the  plasma  part  (which  is  usually  100  to  1000  times  denser  than  the 
beam  part)  requires  a lot  of  simulation  particles;  otherwise  the  noise  from 
the  plasma  is  high  enough  to  obscure  the  instability  itself.  One  example  of 
this  crucial  effect  can  be  observed  in  the  simulation  of  the  so-called  velo- 
city-space ring-plasma  (flute-like)  instability,  which  arises  due  to  the 
interactions  of  energetic  beam  ions  (a  ring  in  velocity-space)  with  cooler 
Maxwellian  plasma  ions,  carried  out  through  a particle  code  by  C.  K.  Birdsall 
et  al.1  Upon  the  suggestion  of  A.  B.  Langdon,  we  constructed  a hybrid  code 
(the  idea  itself  is  by  no  means  new)  using  a particle  code  (ESI  by  A.  B. 
Langdon)  for  the  simulation  of  beam  particles  and  a fluid  equations  for  the 
representation  of  the  plasma  component.  One  simple  example  of  a fluid  Is 
an  Eulerian  fluid,  linearized  (EFL).  The  linearization  is  to  be  monitored 
but  generally  it  is  justifiable  since  the  plasma  component  stays  linear  in 
many  beam-plasma  type  interactions.  Another  possible  fluid  would  be  a non- 
linear Lagrangian  fluid  (LAF)  suggested  by  W.  M.  Nevins.  These  two  fluid 
schemes  will  be  discussed  together  with  successful  applications  of  a hybrid 
code  combining  ESI  and  EFL  (two  versions:  electrostatic  and  electromagnetic) 
to  the  study  of  linear  and  nonlinear  evolution  of  the  above  ring-plasma 
instability.  The  use  of  a hybrid  for  a beam-plasma  type  instability  has 
merits  of  eliminating  the  plasma  particles  as  well  as  the  noise  from  them. 

For  computing  time  on  a CDC  7600  (LLL  MFECC  "A"  machine)  the  hybrid  (ESI  + 

EFL)  takes  about  O.lA  sec  per  time  step  (*  16  usee  per  time  step  per  particle) 
and  this  could  be  further  improved  through  optimization. 

I.  STRUCTURE  OF  HYBRID  CODES 

The  electrostatic  hybrid  (ESI  + EFL)  consists  of  ESI  for  the 
beam  and  EFL  for  the  plasma.  For  EFL.  the  linearized  two-fluid 
equations  in  finite-difference  form  is: 
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The  leap-frog  scheme  is  time  and  space  centered  as  in  Fig.  1(a).  The 
electromagnet ic  hybrid  code  (EMI  +EFL)  is  constructed  by  combining 
EMI  (relativistic  and  nonrelat ivistic  e lect romagnet i c particle  code7) 
with  the  properly  modified  e I ec t romagnet i c version  of  EFL , whose  leap-frog 
scheme  shown  in  Fig.  1(b)  is  similar  to  the  electrostatic  EFL  iFig.  1(a)] 
except  for  some  additional  quantities  kept. 
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Fig.  1(a) 


APPLICATION  OF  A HYBRID  CODE 


,"'1 — I r 


Fig.  1(b) 


The  use  of  ESI  + EFL  in  simulation  of  an  unmagnetized  beam-plasma 
instability  yields  results  (linear  and  nonlinear)  almost  identical  with 
a particle  code.  A more  important  application  of  this  code  is  to  the 
above  ring-plasma  instability.  Simulations  using  a hybrid  code  (marked  x 
in  Fig.  2)  provide  remarkable  verification  of  the  linear  Vlasov  theory  (con- 
tinuous curves  in  Fig.  2)  especially  when  compared  to  particle  code  results 
(Ref.  1,  marked  0 in  Fig.  2).  This  excellent  agreement  extends  to  a 
wide  range  of  plasma  parameters  whether  electrostatic  or  electromagnetic. 

Hybrid  simulations  also  provide  extensive  information  about  non- 
linear evolution  of  this  instability  such  as  beam  spreading,  average  slowing 
down,  and  saturation  level,  which  agree  nicely  with  some  analytic  explana- 
t i ons . 
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III.  LIMITS  ON  THE  LINEARITY  OF  THE  FLUID  PLASMA 


The  linearity  assumption  is  good  as  long  as  the  perturbed  plasma 
velocities  (vi ) remain  small  compared  to  the  wave  phase  velocity  (m/k) . 

Fig.  3 shows  that  this  criterion  is  satisfied  in  many  cases  except  for  strong 
beam  cases  (e.g.,  R = nbean/n plasma  s 1 to  ’)  near  saturation  times.  In 
order  to  follow  the  nonlinearity  of  the  fluid  plasma,  changing  to  the 
Lagrangian  fluid  (LAF)  model  will  be  discussed. 
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Fig.  2 


*C.  K.  Bi rdsal 1 

p.  111*8,  1976. 

2A.  B.  Langdon, 
M.  A.  Mo strom, 
18,  1*70,  1975. 


M.  J.  Gerver  and  N.  Maron,  Bull.  Am, Phys.  Soc.  2 1 , 

Ph.D.  Thesis,  Princeton  University  (1970);  B.  I.  Cohen, 

0.  R.  Nicholson,  A.  N.  Kaufman  and  C.  E.  Max,  Phys.  Fluids 


A paper  will  be  prepared  and  submitted  to  the  Journal  of 
Computational  Physics  next  quarter,  which  will  consist  of  the  detailed 
code  algorithm  with  two  major  applications,  namely  to  the  magnetized 
ring-plasma  instability  and  to  the  unmagnetized  beam-plasma  instability 
[specially  the  off-peak  saturation  effects  presented  in  the  QPR  II 
(July,  1977),  pp.  5-1^) 1 - 
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E. 


lid  PARTICLE- FLU  ID  HYBRID  MODEL  FOR  THE  DRIFT-WAVE  INSTABILITIES 
Yu-Jiuan  Chen 


A simulation  method  is  proposed  for  the  study  of  the  lower- 
hybrid-drift  instability  and  the  mirror-drift-cone  instability  in  the 
sheath  region.  For  simplicity,  it  is  assumed  that  Tg«  T. , and  the 
finite  plasma  beta  effects  (electromagnetic  effects  and  VBq  electron 
orbit  modification)  are  negligible.  Moreover,  primary  emphasis  is 
placed  on  high-frequency  perturbations  characterized  by 


0)  . ■«  | a)  + i y | 3j  (a>  a).) 

ci  1 r 1 ^ ce  c i 
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= 0)  , « 0) 

h ce 


-2  2 -2 

a . « k s a 
i i e 


2 2 
a « L 

e 


where  io  = |e  |B  /m  c is  the  gy rof requency , a = (2T  /m  u )2  is  the 
cs  1 s 1 o s _]  3n  s s s cs 

local  Larmor  radius,  and  L =—  I /n  is  the  density  scale  length. 

3y  Yo 

Therefore,  we  can  employ  the  local  approximation  in  the  density  gradient 
direction  and  treat  electrons  as  a strongly  magnetized  fluid  and  ions 
as  unmagnetized  particles. 

In  general,  the  initial  plasma  configuration  is  illustrated 
in  the  slab  geometry  as  shown  in  Fig.  1.  We  consider  perturbations 
localized  about  y=yQ,  and  the  local  equilibrium  magnetic  field  primar- 
ily in  the  z-direction  with 

?(?)  = e B (y J + e B (yJ  ~ i i-2- 


in  the  vicinity  of  y = y , where  L = ( 3 B / 3y ) | /B  (y  ) is  the  length 

o s x Vb  o o 

scale  characterizing  the  magnetic  shear.  The  electron  drift  across  the 
magnetic  field  B^^  with  the  mean  fluid  velocity 

v = v_  +v . , (**) 

xe  E de 

where  v_  = cE°  / B is  the  |*B  drift  velocity,  v.  = (T  /m  o>  fw(nT  ) | 

t y o de  e e ce  oy  e y^ 

is  the  electron  diamagnetic  drift  veloci  ty,  and  n (y)=n£ (y)  ^ n . (y)  is  the 
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density.  For  the  electron  fluid  drifting  through  a stationary  ion  back- 
ground, the  electron  E*B  drift  motion  and  the  inhomogeneities  in  n, 

B,  and  T provide  the  free  energy  to  drive  instabilities. 

As  L -*■“»,  i.e.,  no  magnetic  shear,  a lid  electrostatic  particle 
fluid  model  (ESI  + FLUID)  can  be  used  to  simulate  these  two  kinds  of 

instabilities  in  E*B  direction,  i.e.,  x direction,  instead  of  2id 

1 ° 

particle  simulation.  A lid  magnetized  electron  fluid  drifts  with 

velocity  v through  a local  Maxwellian  Id  ion  background  for  the 
*e 

lower-hybrid-drift  instability,  or  loss  cone  Id  ions  for  the  mirror- 

drift-cone  instability.  The  electron  E * B drift  and  the  inhomogeneities 

are  taken  care  through  the  drift  velocity  v . Since  there  is  no  boun- 

*e 

dary  effect  in  the  x direction  for  the  slab  geometry,  the  periodic  boun- 
dary condition  usually  used  in  ESI  code  is  appropriate  for  this  case 
(shown  in  Fig.  2).  If  the  magnetic  shear  is  considered,  this  model  can 

be  generalized  to  2£d  (x,  y,  v , v , and  v ). 

x y z 
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Fig.  1 Slab  Geometry  and  Coordinate  System 


1D.  Winske  and  P.  C.  Liewer,  "Particle  Simulation  Studies  of  the  Lower 
Hybrid  Drift  Instability",  Phys.  Fluids  21(6) , pp.  1017-1025  (1978). 


periodic  boundaries 


Fig.  2 lid  Particle-Fluid  Hybrid  Model 

r.i.-  * ■ 

F.  RJET,  PDP  11/04 
Dr.  W.  M.  Fawley 

With  the  exception  of  the  Versatec  controller  unit  described 
below,  all  necessary  hardware  has  been  installed  in  the  UCB  remote  user 
system.  The  PDP  11/0**  has  been  moved  and  secured  into  place  away  from 
the  terminal  rooms  119ME  and  1 1 9MD  where  its  fan  noise  would  have  inter- 
fered with  mental  concentration.  The  DZ-11A  multiplexer  has  been  instal- 
led and  the  8 terminal  lines  run  to  the  proper  rooms. 

In  order  to  run  our  Versatec  electrostatic  printer  off  the 
PDP  11/0**  directly,  it  was  necessary  to  order  in  mid-June  a PDP  ll/DMA 
controller  Model  #121  from  Versatec.  Funds  for  this  acquisition  were 
requested  from  DOE  in  April  and  received  in  June  1978.  An  extra  long 
60 ' cable  was  also  ordered  to  run  from  the  PDP  11/04  to  the  Versatec  unit. 
Delivery  is  expected  by  late  August. 

The  software  necessary  to  run  the  PDP  11/04  as  a terminal  con- 
centrator still  has  not  been  delivered  from  Livermore.  The  primary  diffi- 
culty appears  to  be  the  DEC's  DMC11  network  links  not  working  according 
to  specification.  Consultation  between  LLL  and  DEC  is  in  progress;  until 
the  software  arrives,  the  UCB  PDP  11/04  must  sit  idle. 
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We  have  initiated  efforts  to  obtain  an  LBL  extension  phone 
so  that  we  may  use  their  tie- line  to  LLL.  When  the  UC  system  tie- line 
to  LLL  disappeared  in  February,  our  phone  bill  shot  up  to  >$IOOO/month 
An  LBL  extension  will  cost  <$20/month  and  permit  constant  use  of  the 
LLL  computer.  We  have  received  authorization  from  DOE/San  Francisco 
for  the  extension  and  hope  to  obtain  it  by  mid-August. 

Station  11/ 

PLASMA  SIMULATION  TEXT 


Revisions  of  Part  I have  been  made,  to  be  completed  in  Sep- 
tember. Part  II  revisions  are  progressing  well. 


Station  V 

SUMMARY  OF  REPORTS,  TALKS  AND  PUBLICATONS 
IN  THE  PAST  QUARTER 

Four  papers  were  presented  at  the  Eighth  ConfitA-tnat  on  NlumaZcclZ 
SimuZation  0($  ?lcU> ma6  at  Monterey,  June  28  - 30  (the  PAoattdingA  published 
up  to  four  page  abstracts),  as  follows: 


(1)  J.  K.  Lee  and  C.  K.  Birdsall,  "Particle-Fluid  Hybrid 
Codes  Applied  to  Beam-Plasma,  Ring-Plasma  Instabilities". 

(2)  A.  Peiravi  and  C.  K.  Birdsall,  "Self-Heating  of  Id 
Thermal  Plasmas:  Comparison  of  Weightings;  Optimal  Parameter 
Choices". 

(3)  W.  M.  Nevins,  "Diagnostics  — The  Business  End  of 
a Plasma  Simulation  Code". 

(A)  Y.  Matsuda,  W.  M.  Nevins  and  M.  Gerver,  "A  New  Particle 
Simulation  Model  for  a Bounded  Plasma". 

Prof.  C.  K.  Birdsall  chaired  one  of  the  six  oral  sessions. 
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